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Introduction 

A class of one-dimensional premixed flames that has the potential of eliminating heat loss and stretch effects is the 
burner-generated cylindrical and spherical flames, where a premixed gas flows through a porous cylindrical or 
spherical burner and the resulting flame is formed at some dis tan ce from the burner, as shown in Figs, la and lb. To 
demonstrate that such flames are curved but not stretched, we consider the general expression for the stretch rate 
experienced by a flame surface [1*2], 


=- n -V x(u f x n) + (s f -n)(V-n) ( 1 ) 

where A is the area of an infinitesimal element of the flame surface, u f the flow velocity evaluated at the flame, % 
the velocity of the flame, and n the unit normal vector of the surface. Equation (1) shows that, since the flame is 
stationary (Sf = 0), and the flow is normal to the flame surface (u f x n = 0), the stretch rate vanishes identically. 
Furthermore, while such a flame can be stabilized by heat loss to the burner, the divergent nature of the flow 
provides an additional stabilization mechanism. That is, as the flow rate is continuously increased, the flame recedes 
from the burner and heat loss can be eventually reduced to become negligible as compared to the heat generated at 
the flame. At this point the flame temperature is expected to be close to the adiabatic flame temperature. Further 
increase in the flow rate is then accommodated by an increase in the flame standoff distance to a position where the 
flame speed again balances the flow velocity. Thus, the cylindrical/spherical geometry offers a unique advantage for 
the study of adiabatic, one-dimensional, stretch-free flames that is not present in other commonly used experimental 
systems. Such flames, however, can only be established, and studied, in a microgravity (Jig) environment because of 
the distorting nature of the buoyant flow under normal gravity conditions. 

In view of the above considerations, an experimental and theoretical program chi cylindrical and spherical 
premixed flames in microgravity has been initiated. We are especially interested in: (1) assessing heat loss versus 
flow divergence as the dominant stabilization mechanism; (2) understanding the effects of flame curvature on the 
burning intensity; and (3) determining the laminar burning velocity by using this configuration. In the present study 
we have performed analytical, computational, and pg-experimental investigations of the cylindrical flame. The 
results are presented in the following sections. 

Asymptotic Analysis and Qualitative Behavior 

The problem analyzed is shown in Fig. la. Here a premixture of temperature T u , mass flow rate m, and consisting of 
a deficient, rate-limiting reactant F (fuel) with mass fraction Y F u issues from a unit-length, internally-cooled porous 
cylinder which has a radius r, and surface temperature T„ where the subscript s designates the burner surface. From 
continuity, the mass flow rate m=27trpu is a conserved quantity along the radial distance r, where u is the r-velocity. 
Upon ignition, a cylindrical flame is established concentric with the burner, with the thin reaction zone located at r f 
and the preheat zone extending from the burner surface to r f , where the subscript f designates the reaction zone. 
From the experimental point of view, the parameter which can be independently and conveniently controlled is m. 
Thus for a given m, the cooling rate of the burner and hence the heat loss rate from the flame, Q,, can be adjusted to 
achieve a constant T,. The flame temperature T f and flame standoff location r f will also adjust accordingly for the 
given m and T,. 
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A straightforward application of activation energy asymptotics yields the following flame responses, 

(tr- 


T f =(Tf° + T,-T tt )-q s 


- — m 

<1*=% 


where A(T f ) = 2Da (T^/t,) 2 exp(-t t /T f ) is a reduced DamkOhler number. 


1 l Wo A Vc p 


is the DamkOhler number, Le the Lewis number, m=m/(2nX/c^) , T= c^T/Yp ,,^ Y F =Y F /Yp >u , p= p/Pb. f =r/r 5 , 
q s = [r s /( Y F}1 q c )] q s , T, is the activation temperature, p the density, Y 0 tne mass fraction of the abundant species, X 
the thermal conductivity, c p the specific heat, ^ the chemical heat release per unit mass of fuel, (j, the specific heat 
loss per unit length of burner, and B the frequency factor. The subscripts u and b respectively designate the unbumt 
and burnt states of the flame, and the superscript “o” designates the adiabatic, one-dimensional planar flame. We 
shall also assume that there is no downstream heat loss such that T f = T b . 

For large values of f f the above solution approaches the adiabatic limit, with q,— >0, Tf-iT® and 
A( T f ) -> A(T f °) which is a constant The last relation, when interpreted in die context of Eq. (2), implies that (m/ij) 
becomes a property of the mixture in this limit being independent of the burning rate m and the flame curvature 
which is proportional to if - *. This result has the following implications. 

First if- in, the flame accommodates an increase in the mass flow rate from the burner by moving to a 
larger radial distance and hence attaining a larger flame surface area. The flame is therefore stabilized through flow 
divergence. This stabilization mechanism and hence flexibility do not exist for the flat-bumer flame, which loses its 
planar nature when the burner discharge rate exceeds that of die adiabatic valu 

Second, if we evaluate mat the reaction zone, then (m/r f ) = (p^J/QJCp) is simply a nondimensional mass flux 
at r f , say f f , where f = bj(k/c^) is the nondimensional f. Furthermore, since this quantity is independent of r f and 
hence the flame curvature in the adiabatic limit, it must be equal to the (constant) mass flux of the adiabatic planar 
flame, f°= f£= fj°, for which r, can be considered as an arbitrary length for nondimensionalization. We therefore 
have the following sequence of identities in the adiabatic limit: 


Ifj =<V =tf°) =Led. w 

The inva riants , of ff to flame radius variations provides a simple relation for an experimental determination of 
the laminar burning velocity of die one-dimensional planar flame. That is, since m (=2jcrjP b Uf) is given while the 
flame radius can also be easily measured photographically, the laminar burning flux is simply 


Furthermore, since f°= fjf= p u i£, the laminar burning velocity is then given by s£ = u£= f°/p u . Our experimental 
result agrees well with Eq. (7), as will be shown later. 

Third, since ff= while f f < f n because of flow divergence, we must have £ < f „. Thus if f u is identified as die 
flame, pro pagation rate relative to the freestream, then the concave nature of the present cylindrical flame renders it 
to have a higher flux and hence propagation rate than the planar flame. This is the same curvature effect which 
causes the tip of a unity Lewis number Bunsen flame to have a higher burning velocity than its shoulder. 

The characteristics of the cylindrical flame can be better demonstrated with the following quantities, 


'it) 
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( 9 ) 




where m s is the ratio of the mass burning rate of the cylindrical flame to that of the adiabatic planar flame with the 
same cylinder surface area, orm°= We also define a flame standoff distance df= r £ - r„ with r, * 0 for the 
planar case, and nondimensionalize it by die laminar flame thickness 8° = (X/cJ/f° such that d f = <l f /8 > . 

For comparison, we note that for the nonadiabadic planar flame stabilized over a flat burner [3, 4, 5], 
m 5 = f= frf° such that 



( 10 ) 


Qj = ihj exp(-m s df) (11) 

while Eq. (3) still holds for die relation between T f and q, 

Calculated results using representative values show that both the planar and cylindrical cases yield a dual 
solution behavior in that there are two values of die mass burning rate m s for either a given heat loss rate Q, or a 
given flame standoff distance df The behavior for the planar case is a well known observation [6, 7, 8] which has 
also been recently explained [4, 5] by noting that by taking m as the independent parameter with T s fixed, the flame 
responses in terms of Q, and df are actually unique. By similar reasoning we then expect that the responses of the 
cylindrical flame are also unique' when m s is die independent parameter. 

For the cylindrical flame, it is further found that when die cylinder radius is large, the flame behavior follows 
closely that of the planar flame up to m, = 1, as expected. In this regime heat loss is the dominant flame 
stabilization mechanism. For m, > 1, the planar flame does not exist while the cylindrical flame_can_ still be 
maintained due to flow divergence. In this regime the flame is practically adiabatic, with Q,= 0 and Tf= f°. For 
smaller cylinder radius, the response of the cylindrical flame deviates significantly from that of the planar flame. 
The flame is almost adiabatic aisd is stabilized by flow divergence at higher values of m g . 

To identify the influence of curvature on die burning flux at the upstream boundary of the flame, we have also 
computed ^ = 4/^u°* with the unbumt state evaluated at the location of the 1% temperature rise from the freestream 
temperature according to the suggestion of Tien and Matalon [9]. The results are quite insensitive to the actual state 
selected as long as it is close to the ambient state. It is found that with increasing m,, ^ initially increases to values 
of fj, > 1. It then decreases and eventually approaches the adiabatic, planar case of 4 = 1 as °°. The important 
point to note is that f u can exceed f D °, in agreement with the previous discussion of curvature-enhanced burning flux 
and hence burning velocity. 


Microgravitv Experimentation 

The experimental apparatus for the cylindrical flame is shown in Fig. 2. The cylindrical burner was made of a 
porous bronze tube with 5 tun pore size, 1.25 cm diameter, and 3.76 cm active length. The thickness of die porous 
wall is 4.7 mm. Since uniformity of the flow through the porous wall depends on the pressure differential across it, 
the reactant mixture was supplied from both ends of the burner tube so as to minimize the internal pressure gradient 
in the axial direction of the cylinder. Furthermore, stainless steel tubes with small holes at various locations woe 
positioned inside die burner to further even die axial pressure distribution. A cross sectional view of the burner with 
the internal tubes is shown in the inset of Fig. 2. The axial pressure distribution was measured (under normal gravity 
conditions) by drilling small holes in the porous wall and connecting water manometers to the pressure taps. For 
the best achieved design of the number and placement of the holes in internal tubes at a typical flow rate, the inlet 
gage pressure to the burner was 60 cm water and the axial pressure variation was +/- 2 cm water, a 10% variation. 

Due to the s mall burner size and the limited duration of the microgravity experiments, no cooling was used for 
the burner. A thermocouple was inserted inside the burner wall to record the burner temperature. A hot wire 
mounted on a rotary solenoid was used to ignite the mixture. Experiments were performed under microgravity 
conditions at the NASA-Lewis 22 second drop tower. Mass flow rates were measured by using sonic flow nozzles 
whose upstream pressures were recorded by using pressure transducers. A video camera was used to record the 
flame shape and standoff distance. 
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Computational Simulation 

The experimental situation was numerically simulated by using the Sandia premixed flame code [10], with 
mnHifjpationg to ensure strict mass conservation in the divergent flow field. The program solves the governing 
w piarinnc of mass, energy and species using finite difference approximation, allowing for detailed chemical kinetics, 
multicomponent di ffusion and variable thermal properties. The chemical mechanism used was that developed by 
Egolfopoulos et al. [11], which has 30 species and 153 reaction steps. 

The independent variable for the burner stabilized flame is the mass flow rate m. The necessary boundary 
conditipng are the burner surface temperature, T„ and the reactant mixture mass flux fractions. The boundary 
conditions imposed on the product side are vanishing gradients for species and temperature, dYj/dr and dT/dr. The 
program solves the governing equations to determine the temperature profile through the flame, from which the 
flame location and heat loss to the burner can be found. 

Experimenta l and Computational Results 

A mixture of CH 4 + 20 2 + 7.52N 2 + 5.665He, with an effective Lewis number of 1.88, was used for the cylindrical 
flame experiments. Helium was added to enhance flamefront stability. There are two causes for the enhanced flame 
stability with helium addition. First, the effective Lewis number of the mixture increases, which reduces the 
thermal-diffusional instability. Second, the adiabatic flame temperature is reduced, which reduces the density jump 
across the flame and the propensity to develop the hydrodynamic (Landau) instability. For this mixture the adiabatic 
flame temperature for the unbumed mixture temperature of 300 K is 1885 K, and the calculated laminar flame speed 
is 20.8 cm/s, which corresponds to a mass burning flux of 0.016 gm/cm 2 -s. 

Figure 3 compares the calculated and measured flame standoff distance as a function of the normalized mass 
flow m s . The inset shows a typical recorded flame configuration, which exhibits some large-scale, weak wrinkling. 
This wrinkling is believed to be caused by the pressure variation in the burner and leads to a variable standoff 
digtanee. of +/- 15% in the axial direction. The spread in the experimental data therefore represents the extent of 
flame, nnnnnif ormity as measured by the mimimnm and maximum standoff distances. The flame location for the 
numerical results is taken as that of tire m a ximum CH radical concentration, since this species is chiefly responsible 
for the luminosity of the flame and assumes its maximum concentration in the active reaction zone of the flame. The 
heat loss to the burner is also presented for the calculated results, and shows that around a normalized flow rate of 
m, - 1.4, the heat loss is essentially zoo. At this point the stabilization mechanism has changed from heat loss to 
flow divergence. The figure also shows that all the microgravity data are in the regime where the flame is stabilized 
by flow divergence. This is further substantiated by the experimental result that the thermocouple temperature 
basically remained at the ambient value during the experiment. The comparison shows that the flame is situated 
further away titan calculated, with about 1 to 2 mm differences. Overall, the experimental and computational results 
agree quite well, at least qualitatively, considering the difficulties usually associated with microgravity 
experimentation. 

Figure 4 compares the calculated and observed mass flux through the reaction zone. As discussed, the 
calculated result in the adiabatic limit should be close to the independently calculated laminar burning flux for the 
aHiahatir p lanar flame. Exact agreement is not expected because while the calculated result is a constant for the 
p lanar flame, it depends on the present choice of the maximum CH radical location as that of the reaction zone for 
the cylindrical flame. The uncertainty, however, is not expected to be large because the flame thickness is much 
gmaiw than its standoff distance. It is seen that the calculated laminar mass flux attains a value of about ff = 0.0156 
gm/cm 2 -s in the divergence-stabilized regime, which compares very favorably with the independently calculated 
value of f° =0.016 gm/cm 2 -s for the adiabatic planar flame. Figure 4 further shows that the experimentally measured 
values also agree well with the computed results, being slightly smaller. This demonstrates the feasibility of 
determining the laminar burning flux, and hence the laminar burning velocity, by using the cylindrical/spherical 
flame as proposed herein. 

The difference between the calculated and measured values for the flame location and burning flux can be 
attributed to the small influence of gas-phase heat loss in the experiment, which is not included in the calculations. 
Perhaps most importantly, while the calculation assumes an adiabatic downstream, in the experiments the 
downstream temperature is that of the ambience. This downstream heat loss reduces die flame burning rate, which 
in turn causes the slightly weakened flame to move to a larger radial location so that stabilization is achieved with 
the reduced radial velocity. 
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Concluding Remarks 

In the present investigation we have demonstrated the usefulness of the stationary cylindrical flame for fundamental 
flame studies, especially when the flame is stabilized by flow divergence and hence is basically adiabatic. Since 
these curved flames are stretchless, the influence of curvature on the flame structure and response can be 
unambiguously studied. Specifically, we have shown that the downstream burning flux of such flames are 
unaffected by die flame curvature and hence is equal to the constant burning flux of the adiabatic, freely-propagating 
planar flame, while the upstream burning flux exceeds it due to flow divergence. This understanding is usefiil to the 
study of the general nature of stretched flames subjected to the influence of flame curvature, especially within the 
context of laminar flamelets in the modeling of turbulent flames. 

We have also successfully produced a nearly-cylindrical and adiabatic flame in the microgravity environment of 
the drop tower. While the quality of the flame still needs to be improved, the viability of the concept is nevertheless 
established, especially in light of the favorable agreement between die microgravity results and detailed numerical 
computation. 

A further contribution of the present study is the possibility of determining the laminar burning velocities of 
combustible mixtures by using the cylindrical flame in the nearly-adiabatic regime. While microgravity experiments 
are fairly complex and possibly also costly in terms of their operation, the simplicity of the experimental 
configuration and concept renders this methodology a useful substantiating alternative to other approaches for such a 
determination, at least for reasonably smooth flames. 

Finally, we do wish to emphasize again the difficulty in obtaining perfectly cylindrical flames without wrinkles, 
even with the addition of helium. If such wrinkles are inherent to the flame, as manifestations of some 
hydrodynamic instabilities [12], then there would exist a limitation on the quantitative accuracy in determining the 
flame responses such as the burning rate of practical combustion mixtures. Studies are currently underway to assess 
the nature of such wrinkles. Experiments will also be conducted with spherical flames; it is anticipated that the 
additional curvature could partially “absorb” the extent of wrinkling. 
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Figure 1. Configurations of (a) cylindrical, and 
(b) spherical flames which are curved but unstretched. 



Figure 3. Computed heat loss flux and flame standoff 
distance, and experimental flame standoff distance, as 
functions of the discharge rate. Inset is a representative 
experimental flame configuration. 



Figure 4. Computed and reaction zone mass flux as 
function of die discharge rate. 


Figure 2. Schematic of the experimental apparatus. 
Inset is the cross-section of the cylindrical burner. 
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